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ТЕХНИЧЕСКИЕ НАУКИ 


УДК 613.6:621.43 
Б.Ч. МЕСХИ, А.Н. СОЛОВЬЕВ, Ю.И. БУЛЫГИН, Л.Н. АЛЕКСЕЕНКО, Е.И. МАСЛОВ 


МОДЕЛЬ РАСПРОСТРАНЕНИЯ ВРЕДНЫХ ВЕЩЕСТВ 
И ИЗБЫТКОВ ЯВНОЙ ТЕПЛОТЫ В ПРОИЗВОДСТВЕННЫХ ПОМЕЩЕНИЯХ 


Предложена нестационарная трехмерная математическая модель, описывающая процессы конвективно- 
диффузионног о тепломассопереноса вредных веществ и избытков теплоты в производственном помещении с 
активной вентиляцией. Модель учитывает как стационарные, так и передвижные источники загрязнения и 
позволяет прогнозировать качество воздуха рабочей зоны в помещениях сложной конфиг} урации. 
Ключевые слова: загазованность, рабочая зона, прогнозирование, распространение вредных веществ, 
моделирование, параметры микроклимата, производственная среда. 


Введение. Повышенная загазованность рабочих зон производственных помещений машино- 
строительных и ремонтно-обслуживающих предприятий сочетается, как правило, с неблагоприят- 
ными параметрами микроклимата (повышенная температура и высокая подвижность воздуха). 
Так, концентрация вредных веществ в рабочих зонах вблизи действующих источников загрязне- 
ния может превышать максимально-разовые ПДК». в десятки и сотни раз. А параметры микро- 
климата в различные периоды года также не соответствуют нормативам. Поэтому не вызывает 
сомнения, что задача снижения загазованности и сохранения здоровья работников таких цехов 
при вредном воздействии токсичных выбросов является актуальной. 

Для решения данной задачи разработана математическая модель, описывающая процессы 

тепломассопереноса вредных веществ в объеме производственного помещения, которая позволя- 
ет определить опасные зоны загазованности и области с неудовлетворительными параметрами 
микроклимата, оказывающие негативное воздействие на работников. 
Математическая модель тепломассопереноса вредных веществ и избытков явной теп- 
лоты. Система уравнений в общем виде. Распространение примесей, движение воздуха 
(в потенциальном приближении) и изменение температуры в рассматриваемом объеме описыва- 
ется следующей системой уравнений [1, 2]: 
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жду коэффициентами турбулентной диффузии помещения; С, - концентрация 1-го вредно- 


где Л, — коэффициент турбулентной теплопроводности, отображающий связь ме- 


го вещества в помещении; и‚у,и -— компоненты вектора скорости воздушной среды; 
у’, — скорость оседания/подъема примеси/газов; А, =(А,,А,,А,) — коэффициенты турбу- 


лентного обмена помещения; О — коэффициенты диффузии; О/ (Е) — интенсивность вы- 


х,у,2 
броса г-го вредного вещества от /-го источника в помещении; 6(’—и,) — дельта-функция 
Дирака; Р -— потенциал скоростей движения воздуха; р — плотность воздуха; Т -— темпера- 
тура; #, ие »Ур»2 а координаты источника выброса; и‚у,м — компоненты скорости; 


С, -удельная изобарная теплоемкость; ^, - коэффициент теплопроводности; /., — коэф- 


фициент радиационной теплопроводности; О/ - интенсивность внутренних источников те- 
пла; Рг, — диффузионное число Прандтля; Рг, —турбулентное число Прандтля. 


Граничные условия для уравнения (1): 
— на боковых стенках 
Я. + 
ди 
В этом равенстве указано, что стенки помещения, в котором проводятся исследования, 
являются непроницаемыми для концентрации загрязняющего вещества; 
— на входной и выходной границе 
0©_ а 


НЫ, 
ди А 


0. 


п 


0%, < 


где а= |. ь. — означает, что поток примеси зависит от ее концентрации, скорости ®, и на- 
т 
правления воздушного потока, индуцированного работой вентиляции. 
Граничные условия для уравнения (2): 
— на твердых стенках 
ОР _ 
би 
где и —единичный вектор внешней нормали. 
Это означает, что стенки исследуемого помещения являются непроницаемыми для воз- 
душных потоков; 
— на выходной границе 


0, 


ОР 
бт тя 
где о, - известное значение скорости; 
— на входной границе (границе втекания воздушного потока в помещение) 
Р = сопз{ (условие Дирихле). 
Граничные условия для уравнения (3) [3]. 
Начальные условия заключаются в задании распределения поля значений температуры в 
начальный момент времени (1=0), т.е. предшествующий расчетному. Они должны быть заданы в 


виде функции Т_, = /(х,у,2) - для пространственной задачи. Граничные условия — тепловые 


условия у поверхности тела, которые задаются в более сложном виде. При решении задач тепло- 
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проводности принято различать четыре часто встречающихся способа задания граничных усло- 
вий, так называемые граничные условия первого, второго, третьего, четвертого и пятого рода. 
Граничные условия третьего рода заключаются в задании температуры поверхности тела и окру- 
жающей его среды и задании теплообмена (коэффициента теплопередачи) между поверхностью 
этого тела и окружающей средой по закону Ньютона. Таким образом, количество теплоты, отда- 
ваемое (или получаемое) единицей поверхности температурой Т, за единицу времени в окру- 


жающую среду температурой Т., прямо пропорционально разности температуры поверхности и 
окружающей среды: а, =а(Т, -Т,) [4]. Количество теплоты, отдаваемое (или получаемое) по- 


верхностью в окружающую среду, должно быть равно количеству теплоты, подводимому к этой 
поверхности за счет теплопроводности, которое определяется по закону Фурье. Приравняв эти 
потоки, получим новое выражение для задания граничных условий третьего рода: 


ОТ а 
—=-— (Г. -Т,). (4) 
дп ^, 
В условии (4) должны быть заданы коэффициент теплоотдачи а и температура окружаю- 
щей тело среды Т. . В ходе экспериментальных работ граничные условия могут быть уточнены. 


При турбулентном течении потока и турбулентном переносе коэффициенты А„ пропор- 
циональны диссипируемой энергии и определяющему размеру рассматриваемого источника. Для 
вентилируемых помещений, зная величину энергии, поступающей в помещение, и определяющий 
размер источника, можно определить коэффициент турбулентного обмена в плоскости (х, у) [5]: 


— ри Узу 4/3 
А, =А, = 0,2551 №, 
где # - кинетическая энергия воздушных потоков в помещении; Г, = определяющий размер по- 


мещения, равный для прямоугольных отверстий и сечений неправильной формы /, =\Е, 


здесь Р’— площадь сечения помещения, перпендикулярная направлению движения воздуш- 

ных потоков. 
Коэффициент турбулентного обмена по высоте помещения находится из зависимости, 
предложенной М.Е. Берляндом для наибольших классов устойчивости воздушной среды 
А, = А, (2/2,), ЧТО допустимо для исследуемого помещения. Здесь А_ — коэффициент турбу- 


лентного обмена на высоте 1 м, 2,=1 м, 2 - текущая высота помещения. 
Кинетическая энергия воздушных потоков в помещении определяется из формулы: 


&=5. +8. +81 —&,с, Куда входят соответственно энергии приточных струй =,„., тепловых 


то 
струй =„., движущихся предметов &и и вытяжных струй &„‹. Величины Ех, 5 
ваются в модели рассеивания вредных веществ и находятся из формул: 
а -. ь 
Я — у_ 2! 


п 


= 20? 1+п 
и Сп ор, 4 

т -Узр т 
2 3600 
где [/у, - кратность воздухообмена в помещении свободным объемом У,; у - средняя скорость 
выхода воздуха из приточных отверстий; & — коэффициент местного сопротивления на вы- 
ходе из воздухораспределителей; о’ - теплонапряженность объема; = - ускорение свобод- 
ного падения; 2 — расстояние от полюса до рассматриваемой точки; и=0,8 — эксперимен- 
тальная константа; С, — теплоемкость воздуха на удалении от источника; 7, — температура 
воздуха в помещении на удалении от источника; р, — плотность воздуха на удалении от ис- 
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точника; & -— коэффициент аэродинамического сопротивления движущегося предмета; 


Е, -— площадь поперечного сечения движущегося предмета; У„ -— скорость движущегося 
предмета; т — средняя продолжительность движения предмета с данной скоростью в тече- 
ние 1ч; С’ - масса воздуха в объеме помещения. 

Так как С’, = Гр, где Г - объем помещения, то 
3 
5 Е И. У п т 
д.П 
27 3600 
Явное теплопоступление в помещение учитывается для источников, оборудованных мест- 
ными отсосами, как 





О, = СьогОог” (Тог” -То А - Фог Ал, 
а для источников ими не оборудованных 
О, = СрогОог" (Гог” — То) Ад 
где стт, с*г — изобарные теплоемкости отработавших газов (ОГ) на режимах максимальной на- 


грузки и холостого хода; Тг“, Тх. - температура ОГ на различных режимах нагрузки; 


А, — количество стендов. 


Для реализации предложенной математической модели разработано программное обеспе- 
чение, реализующее интерфейс, который обеспечивает взаимодействие пользователя с пакетом 
НехРОЕ и базой данных, в которой содержатся все необходимые данные для проведения 
модельных расчетов. Интерфейс системы создан так, что любой пользователь, имеющий элемен- 
тарные навыки работы с М/Ипдом5-приложениями, может без труда работать с программным 
продуктом. 

Результаты модельных расчетов. В предложенной математической модели возможно про- 
ведение расчетов для различных схем организации движения газовоздушных потоков в помеще- 
ниях, а также для штатных и аварийных ситуаций (например, неработающая принудительная 
вентиляция, разрывы отсосов и т.д.). Численные результаты получены на основе конечно- 
элементного моделирования системы уравнений (1)-(3). На рис.1, 2 представлены рассчитанные 
по модели поля относительных концентраций оксидов азота на высоте 2=4 м за один час работы 
обкаточного участка « поля подвижности воздуха в помещении в зависимости от расчетных пе- 
риодов года. 








И =1 Сопоещиавов НЯ уг . р , , 

РР | | . 
и й № =0,05 К =0,74 К =0,21. [ соаавиьвы 
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а) 6) 
Рис.1. Поля относительных значений концентрации МО, в воздухе участка обкатки на высоте 2=4 м в штатном режиме: 
а- «теплый» период года; б-— «холодный» период года 
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а) 6) 
Рис. 2. Поля подвижности воздуха в помещении участка обкатки на высоте 2=4 м в штатном режиме: 
а= «теплый» период года; б-—«холодный» период года 


На рис.3 показана динамика изменения относительных значений максимально-разовых 
концентраций оксидов азота в воздухе помещения на высоте «среза» 4 м в течение одного рабо- 
чего часа. 
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а) 6) 
Рис.3. Динамика изменения относительных концентраций оксидов азота в воздухе помещения 
участка обкатки в течение одного рабочего часа: 


а- «теплый» период года; б- «холодный» период года; а-е - выбранные точки помещения 


Выводы. Показана возможность применения разработанной модели для уточненного 
прогнозирования качества воздуха в рабочей зоне производственных участков; созданное про- 
граммное обеспечение позволяет оперативно определять опасные зоны загазованности, поля 
подвижности воздуха и температур для различных производственных участков; результаты моде- 
лирования могут быть использованы для определения рациональных мест расположения возду- 
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хораспределителей в системах приточно-вытяжной вентиляции помещений и для выбора мест 
размещения датчиков газового контроля. 
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ПЮ ассоипЕ Бой $авопагу апа то Е 5оигсез оЁ ро/ийоп апа а//ои/з {о ргебсЕ ай диаМу оЁ {пе иогКтад агеа тп 
{Пе ргет/5ез о {те сотр/саеа сопйдигаНоп. 

Кеу игогаб: дах сощепЕр иода агеа, Югесазйпд, @зеттпайНоп оЁ Ватт!!! зирапсез, тодета, тисгос/тае 
рагате{егз, иогКта епигоптепЕ 
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МЕСХИ Бесарион Чохоевич (р. 1959), ректор Донского государственного технического универ- 
ситета, заведующий кафедрой «Безопасность жизнедеятельности и защита окружающей среды», 
доктор технических наук (2004), профессор (2006). Окончил Ростовский-на-Дону институт сель- 
скохозяйственного машиностроения (1982). 

Область научных интересов: обеспечение безопасных условий труда; изучение условий труда, 
причин производственного травматизма и профессиональных заболеваний в производственных 
процессах крупных предприятий города и в строительстве; разработка инженерно- 
эргономических требований в области человеко-машинных взаимодействий и их оперативной 
реализации; изучение методов проектирования технических систем безопасности жизнедеятель- 
ности с учетом человеческого фактора. 

Автор более 100 научных публикаций, в том числе 4 монографий. 


гесерНоп@аоп$и.ги 


СОЛОВЬЕВ Аркадий Николаевич (р. 1954), заведующий кафедрой «Сопротивление материа- 
лов» Донского государственного технического университета, доктор физико-математических наук 
(2006), профессор (2007). Окончил Ростовский государственный университет (1976). 

Область научных интересов: механика твердого деформируемого тела, обратные задачи, метод 
конечных элементов, генетические алгоритмы. 

Автор более 200 научных публикаций, в том числе 3 монографий. Имеет 1 патент. 


БУЛЫГИН Юрий Игоревич (р. 1966), профессор кафедры «Безопасность жизнедеятельности и 
защита окружающей среды» Донского государственного технического университета, доктор тех- 
нических наук (2007), профессор (2010). Окончил Ростовский институт инженеров железнодо- 
рожного транспорта (1988). 

Область научных интересов: улучшение условий труда работников транспортных и машинострои- 
тельных предприятий. 

Автор более 100 научных публикаций. 


Буучг_го$юу@тай.ги 


АЛЕКСЕЕНКО Людмила Николаевна, старший преподаватель кафедры «Безопасность жизне- 
деятельности и защита окружающей среды» Донского государственного технического универси- 
тета, кандидат технических наук (2010). Окончила Донской государственный технический универ- 
ситет (2001). 

Область научных интересов: улучшение условий труда работников транспортных и машинострои- 
тельных предприятий. 

Автор более 30 научных публикаций. 


[а1еКзеепко@аоп$и.ги 

МАСЛОВ Евгений Игоревич (р. 1985), аспирант кафедры «Безопасность жизнедеятельности и 
защита окружающей среды» Донского государственного технического университета. Окончил 
Донской государственный технический университет (2007). 

Область научных интересов: охрана труда, информационные системы и технологии. 


Автор 5 научных публикаций. 


Ма)опк@уапаех.ги 
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Везамоп СН. МЕЗКНТ (1959), Кесог оЁ Роп Зае Тесйтса! Упмег®Ку, Неа4 оЁ Пе Ш апа Епм:- 
гоптепё Ргф{есНоп 54епсез Верайтепё, Боп За Тесптса! Упм\ег®Ку. РАБ т баепсе (2004), Рго- 
Геззог (2006). Не дгадиаеа гот Ко$юу-оп-Боп п&йЙЩе ог АдисиКига! Епатееипд (1982). 

Кезеагси пеге${5: есоюду, епутоптепЕ ргоесНоп. 

АцПог оЁ оуег 100 ри сабоп$, тацата 4 топодгарнз. 


АгКа4ду М. ЗОГО\МУЕУ (1954), Неад оЁ 1е 5гепо{Й о Маепа!5 ВерацтепЕ Боп Зае Тесйпгса! 
Упмег®Ку. РАБ т РВу$с$ апа Ма\И$ (2006), РгоГе$зог (2007). Не дгадцаеа гот Ко$юу $а{е Упмег- 
$Ну (1976). 

Кезеагсй т\еге${5: теспап!с$ оЁ зо!а аегогтаЫе Боду, шуегзе ргоМетз, ИпКе ветепЕ тео, де- 
пейс авогЁНт5. 

АцПог оЁ оуег 200 ри сабоп$, паиата 3 топодгарй$ апа 1 раег. 


Уигу Г. ВУЕУСТМ (1966), РгоРе5ог оё {Пе Ше апа Епитоптепе Ргоф{есНоп 54епсе$ Оерайтет\, Боп 
Зае Тесйтса! УпмегзКу. РАО т 5аепсе (2007), Ргогезог (2010). Не дгадиаеа гот Козфом Тп5 ще 
ОЕ Кай Тгапзро": Епдтеейпд (1988). 

Везеагси т\еге{&5: ипргоуетепЕ оЁ мюгКта сопа®оп$ оЁ {гапзро[ ап епдтеейпад могКегс. 

Ац ог оЁ оуег 100 рисабопз. 


ЕчатНа М. АЕЕК$ЕУЕМКО, 5епог [ескигег оЁ {Пе Ше апа ЕпмутгоптепЕ Рг{есНоп 5аепсез ОБерап- 
тепе, Воп З{ае Теснтса! Упмег®йу. Сап4ае оЁ 54епсе т Епатеейпоа (2010). Пе дгадиаеа Нот 
Боп За{е Тесйтса! Упм№ег®Ку (2001). 

ВезеагсИ т\еге{5: ипргометепЕ оЁ мюгКта сопаоп$ оЁ {гапзро[ ап епдтеейпад могКегс. 

Ацог оЕ оуег 30 ри сайопз. 


Еудепу Г. МАЗЬОУ (1985), Ро$югадиае зшиаепе ог {Пе Ше апа ЕпугоптепЕ Рг{есЧНоп Берайтепе, 
Боп За{е Тесйтса! Упм№егоКу. Не дгадцаеа от Боп ${ае Тесйпса! УпмегеКу (2007). 

Кезеагси т\егез(: |абоуг рго{есНоп, шюгтаНоп зу$етз апа {есппооЧе$. 

Ацог оЁ 5 ри Гсайопсз. 
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